      subroutine Flux1_2(ie,F1_2)
	  use global
      Implicit none
	  integer :: ie, kloc_e(2)
	  Real(8) :: L_e,vcor_e(2),Flux_e(2,3),U1_2(2,3),F1_2(2,3),U_Temp(2,3)
	  kloc_e=[ie,ie+1]
	  U1_2=0.D0
	  vcor_e=vcor(kloc_e)
	  L_e=Abs(vcor_e(2)-vcor_e(1))
	  
	  Flux_e=F(kloc_e,:)
	  U_Temp(1,:)=-0.5D0*Delta_t*(Flux_e(2,:)-Flux_e(1,:))/L_e
	  U_Temp(2,:)=U_Temp(1,:)
	  
	  U1_2=U(kloc_e,:)+U_Temp
	  
	  !open(10,file='Test.txt',access='APPEND')
	  !  write(10,*) U1_2(1,:)
	  !close(10)
	  
	  Call Flux(U1_2,F1_2,2)



      end subroutine